Interlayer Quasiparticle Transport in the Vortex State of Josephson Coupled 

Superconductors 
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We calculate the dependence of the interlayer quasiparticle conductivity, aq, in a Josephson coupled 
d-wave superconductor on the magnetic field B || c and the temperature T. We consider a clean 
superconductor with resonant impurity scattering and a dominant coherent interlayer tunneling. 
When pancake vortices in adjacent layers are weakly correlated at low T the conductivity increases 
sharply with B before reaching an extended region of slow linear growth, while at high T it initially 
decreases and then reaches the same linear regime. For correlated pancakes aq increases much more 
strongly with the applied field. 



Experimental study of quasiparticle properties in high- 
temperature superconductors so far has focussed almost 
exclusively on thermodynamic and in-plane transport 
properties. Two remarkable theoretical results connect 
these properties with the symmetry of the superconduct- 
ing gap: a) a massless Dirac spectrum of nodal quasipar- 
ticles in a d-wave superconductor leads to a finite den- 
sity of states (DOS) at the Fermi level in presence of 
impurity scattering [y and to a universal (impurity inde- 
pendent) low-temperature limit of the in-plane thermal 
conductivity, kqq = limT^o '*a6(0,T)/r, in the absence 
of a magnetic field. |g,§l; b) in the vortex state the DOS 
is locally enhanced due to the effect of the supercurrents 
around the vortex cores on the near-nodal quasiparticles 
(Volovik effect) [Q. Consequently, the effect of the mag- 
netic field on the electrical and the thermal conductivi- 
ties is twofold: the net enhancement of the DOS tends 
to increase the conductivity, while local changes in the 
phase and amplitude of the order parameter lead to its 
decrease. ^,|| 

The enhancement of the DOS with the magnetic field, 
B, has been invoked to explain the increase of the spe- 
cific heat with B observed in YBa2Cu307_5 (YBCO). 
M The in-plane thermal conductivity, Kab, has been 
found to have a highly non-trivial dependence on the 
applied field and the temperature in both YBCO and 
BizSrzCaCuaOg+s (Bi-2212) compounds. |J|l At low 
temperatures Kab increases with B, while at T > 5 — 10 
K it decreases and then often saturates in fields above a 
few Tesla. In some zero-field cooled Bi-2212 samples a 
hysteretic field dependence with two plateaus at differ- 
ent values of Kab{B) has been observed. |§] This indicates 
that the behavior of Kab{B) is sensitive to the vortex ar- 
rangement in the sample. 

Franz has appealed to the effects of disorder in 
the vortex positions to explain the observed high-field 
plateau in Kab{B). He has argued that randomly posi- 
tioned vortices act similarly to impurities leading to the 



enhancement of both the DOS and the electron scatter- 
ing. These two effects compensate each other at high 
fields (in analogy to the zero-field result being insensi- 
tive to impurity concentration) leading to the universal 
Kab{B, T)/T -^ Koo- However, to obtain this result Franz 
has made questionable approximations: he has replaced 
the spatial average of the product of two Green's func- 
tions by the product of the averages, and has assumed 
a Lorentzian shape of the distribution function of the 
in-plane supervelocity, which differs from a realistic dis- 
tribution in the vortex state, see below. Clearly, this 
approach cannot explain the different (non-universal) 
plateau values. [|| 

Recently it has been demonstrated |l^,|ll|l that the in- 
terlayer quasiparticle electrical transport can be stud- 
ied directly in the resistive state of small area samples 
(mesas) fabricated from the Josephson coupled super- 
conductors such as Bi-2212. In the energy and temper- 
ature range below « 3 meV, the experimental data are 
well understood in the framework of a Fermi-liquid model 
for the near-nodal quasiparticles in a d-wave supercon- 
ductor assuming (i) clean limit, (ii) resonant (unitarity) 
impurity scattering, and (iii) dominant contribution to 
the c-axis conductivity from coherent (conserving the in- 
plane momentum) interlayer tunneling. pi| ] The authors 
of Ref. [ni| introduced the universal quasiparticle c-axis 
conductivity <jq(T = 0,B = Q). This quantity is insensi- 
tive to the impurity vertex corrections which modify the 
in-plane conductivity as discussed by Durst and Lee. 0] 

So far the experimental data on the magnetic field 
dependence of interlayer quasiparticle conductivity have 
been quite scarce. Measurements of (Jc{B) in high fields, 
where Uc ~ erg, reveal a linear increase of ac with B on 
the scale of 40 T. |l^ For mesas Yurgens et al. [Q have 
reported a weak field dependence for B < 7 T without 
quantitative results for its functional form. 

In this Letter we address the effect of vortices on the 
c-axis quasiparticle conductivity, aq{B,T), as a function 



of the magnetic field B || c, in Josephson coupled su- 
perconductors in the framework of the approach used in 
Ref. [|l|. We find that at low (high) temperatures for c- 
axis uncorrelated vortices aq{B,T) increases (decreases) 
with B, before reaching an extended region of slow linear- 
iii-B increase. In this high field regime the increase in the 
DOS is largely compensated by the enhancement of the 
electron scattering at tunneling, due to disorder in the 
positions of vortices in neighboring layers. 

To calculate the effect of vortices on the interlayer 
transport in the framework of the d-wave model we 
employ a semiclassical approach. M The supercurrents 
around vortex cores lead to a Doppler shift, e„(k, r) = 
k ■ Vs„(r), in the quasiparticle spectrum, £(k, v^) = 
Ek + £«• Here k is the quasiparticle momentum and 
Vs„(r) is the supervelocity at point r inside layer n, 
Ek ~ Kk+^^C*)]^^^ is the quasiparticle energy at B = 0, 
Ck = VF • (k - kp) and A(k) = Ao{kl - kl)/P. As 
c-axis currents are parallel to the field (in contrast to 
the in-plane currents), g the net interlayer transport 
is determined by the spatial average of local transport 
coefficients. In the following we neglect correlations be- 
tween the positions of impurities and those of vortices, 
and average over the distribution of each independently, 
see below. We also neglect the temperature dependence 
of Aq and small effects due to Zeeman splitting. Then 
the interlayer quasiparticle conductivity is given by pH 
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where aq{0,0) = 2e'^t^N{{))srj/TTfilS.Q is the universal c- 
axis conductivity, A^(0) is the 2D DOS, t{k) is the inter- 
layer transfer integral, io — ^(kg), kg are the positions of 
the gap nodes on the Fermi surface, and < 77 < 1 is the 
weight for the coherent tunneling. Further, (...) denotes 
the spatial average, and 

^^ L{cj,Ei,) 

El., 



A{k,^)^ 1 + ^ 



Ek) 2 



is the spectral density averaged over impurities with 
L{uj,E) = j{Lu)/Tr[{E - n{uj))^ +7^(w)]. The functions 
uj — n(uj) and —7(0;) are the real and the imaginary part 
of the self-energy respectively. B In the unitarity limit 
the effective scattering rate of quasiparticles is 7(0;) w 
70 — (jj^/870 when a; ^ 70 and 7(0;) w TrjQ/2\uj\ when 
^ 3> 7o, where 70 ~ {TivoAq)^''^ and vq is the bare scat- 
tering rate. The renormalized frequency is Q{lo) w uj/2 
when w <C 7o and fl{uj) « tj if w ;» 70 . 

At low temperatures T <C Aq the quasiparticle current 
comes mainly from the regions near the gap nodes. In 
the vicinity of a node we can linearize the quasiparticle 
spectrum and obtain from Eq. (|l|) 
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(L(w + e„+i,E)[L{uj + en,E) + il/3)L{Lu + e„, -E)]), 

where e„(r) = kg • Vs„(r) and wc have set t{k) ~ Iq. 
The characteristic energy scale for the Doppler shift, 
e„, is es = hvp/a, where a — (<i>o/-B)^/^ is typi- 
cal intervortex distance. The semiclassical approach is 
valid for es ^ Aq, i.e. for fields B ^ Ba, where 
Ba — *&oAo/?i^w|-. This approach also does not account 
correctly for the quasiparticles in the regions near the 
vortex cores, which leads to corrections of the order of 
Saq/aq « e^/Ag = B/Ba- Another important field scale 
can be obtained by comparing eg with the scattering rate 
7o; the field By — ^Qj^/h'^Vp separates the regimes of 
impurity dominated and field dominated behavior. Using 
the parameters Aq ~ 25 meV, 70 « 2 — 3 meV and 
VF ~ 1.5 - 2.5 • 10^ cm/s [|3| we obtain B-y « 0.2 - 0.6 
T and Ba « 40 - 80 T. 

We now rewrite the spatial average in Eq.(||) as the 
average over the probability distribution of the Doppler 
shift, e, which is fully determined by the vortex arrange- 
ment. In Bi-2212 crystals the 3D vortex lattice is de- 
stroyed by pinning at the "second peak field" B « 0.02- 
0.05 T (see, e.g., [Q), and at higher fields pancakes in 
neighboring layers are only weakly correlated. Conse- 
quently, we consider the limits of c-axis-correlated and 
uncorrelated pancakes. In both cases the average de- 
pends only on the probability distribution function of 
e in a single layer, 'P(e), which is related to the dis- 
tribution P{px) of a single component of supermomen- 
tum px ~ 2mVsx, where m is the effective mass, by 
Vie) = {2/vf)P{Px = 2e/vF)- The function P{px) is de- 
termined by the pancake configuration {R^ = (Xi,Yi)}, 
and, when a is smaller than the London penetration 
depth Aab, is given by P(p2,) = {S{px-lT'Y.i^t/Ri)){Ri}- 
For an isolated vortex p ex h/R where ^ab <^ i? <C a is 
the distance from the pancake center, and ^ab is the co- 
herence length. Consequently, V{e) has a universal tail 
■p(e) = 7re^/(8e^) at es < e < Aq, while its behavior 
at e ^ es depends on the actual positions of vortices. 
However, in absence of correlation between the positions 
of vortices and impurities V{e) involves a single energy 
scale es for any e, and depends on e only via e/es, i.e. 
'P{e)de = 'PB{^/^B)de/eB- In the clean limit we can ex- 
tend the asymptotic behavior 'P(e) oc 1/e^ to infinity 
when the integral over e in Eq. (g) converges. Therefore, 
importantly, up to terms of the order of 70/Ao, the con- 
ductivity depends only on the dimcnsionless parameters 
es/70 = \fBjB^ and T/70. 

For c-axis-correlated vortices e„ « Cn+i, and the aver- 
age in Eq. (||) means (:r(e„, e„+i)) « / dtP(t)T(e, e). At 
low temperatures, T <gc 70, both in weak (B <^ _B-y) and 
in strong (B ^ B^) fields the leading field-dependent 
part of the conductivity varies as (e/70)^. We therefore 
separate aq into the contributions from the regions e < 70 
and e > 70, and use the asymptotic behavior of ^(u) and 
7(0;) in each case to estimate 



a,(i?,0)-a,(0,0) ^ l (e")e<^o Stt^ (e^) 
a, (0,0) ~6 7o' 8 7o' 



g>70 



(3) 



As the average is dominated by the tail of 7'(e) we can 
evaluate aq in weak and strong fields with logarithmic ac- 
curacy. Moreover, we interpolate between the two limits 
to obtain also a qualitative description at B ^ B. 
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Therefore at low T for the 3D-ordered case aq increases 
quasi-linearly with B in the entire field range up to B^ ■ 
For the c-axis uncorrelated vortices the average in 
Eq. (0) has to be taken independently in each layer, 
(J'^(e„,e„+i)) = j d€iV{ei) j de2'P{e2)T{ei,€2). Again, 
for T <C 7o and B <^ Bj we expand in en/70: a-nd obtain 
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At high fields the variation of e„ is of the order of eb ^ 70 
and consequently {L{uj + e„, E)) « P{E — uj). Then for 
B^ -^ B <Si Ba the conductivity is given by 



<Jq{0,0) 



B 1"°° 

Ci + — , Ci = 2^2 / deeV^ie), 
Bq Jq 



(6) 



The most important contribution to the linear, in B, 
term is due to the increase in the tunneling away from 
the nodes, i(k) ~ ip +ii0^, where (p is the angle between 
k and the nodal direction, and ii ^ to, ll^ which yields 
Bo ^ {to/ti)BA- Other contributions, due to deviations 
of the quasiparticle spectrum from the massless Dirac 
form of the linearized dispersion and due to corrections 
to the semiclassical approximation in the vicinity of the 
vortex cores, only enhance cr, on the scale of Ba S> Bq. 
Very importantly, due to scaling of V{e) = 'Pb{^/^b)/^b 
(see above), Ci in Eq. (^) is a constant which depends 
solely on the shape of the distribution 'P(e), but not on 
the magnetic field. If 'Pb{x) is monotonous Ci depends 
on its asymptotic decay. For the Lorentzian (when vor- 
tices act exactly as the impurity scatterers do at S = 0) 
Vb{x) oc a;~2 at x 3> 1, and Ci = Cl = 1, while 
for the Gaussian distribution Ci = Cq = n/2. In the 
vortex state Vsix) oc a;~^ at large x, and we expect 
1 < Cl < 7r/2. Therefore to find aq we need accurate 
information on pancake arrangement to determine the 
supervclocity distribution 'P{e). 

Below the irreversibility line Tirr{B) at high fields and 
in the vortex liquid state pancakes do not possess long 
range order inside layers due to pinning and thermal fluc- 
tuations. To calculate aq{B,T) in these regimes we use 
the distribution function P{px) obtained by numerical 
simulation of the 2D pancake liquid at different values 
of the dimensionless temperature t = T^ff /2'kEo, which 



characterizes pancake disorder inside the layers. Here 
Bo = ^qs/Wtt^X^i^ is the characteristic energy of pan- 
cake interaction. In the liquid state T^ff = T, while in 
the glass state it is reasonable to take Te// ~ Tirr{B). 
The calculated function P{px) is shown in Fig. 1 for sev- 
eral values of t; the inset shows that the parameter Ci 
grows with t from 1.08 to 1.22. The full field dependence 
of aq{B,T)/aq{0,0) at T <C 70 obtained from these dis- 
tribution functions is shown in Fig. 2. In computing Og 
we have determined 7(w) and fl{uj) self-consistently llO] 
in the unitarity limit, choosing the scattering rate so that 
70 = O.lAo, typical of Bi-2212 samples. |l| 

Due to short range correlations inherent to strongly in- 
teracting vortices at i3 > $o/A^^ the calculated function 
P{px) differs significantly from the Gaussian (obtained in 
Ref. |l7| assuming uncorrelated and random vortex po- 
sitions). The distribution function depends on a single 
variable e/e_B only if the positions of vortices and those of 
impurities are uncorrelated; this is the case for a vortex 
solid and for the liquid phase in presence of weak pin- 
ning. In general, below the irreversibility line impurities 
act as pinning centers, and the distribution function de- 
pends on two variables, e/e^ and B/^oni, where rii is 
the impurity concentration in a layer. We believe nev- 
ertheless that in that regime our approach gives at least 
the correct qualitative behavior. 

At T ;3> 7o, we also consider the intermediate fields, 
B^ ^ B ^ Bt — ^oT'^/fi^Vp, when the quasiparticle 
concentration is determined by the temperature and the 
primary effect of vortices is to increase the scattering at 
tunneling. Consequently aq decreases with field 



(rq{B,T) 
(7g(0,0) 
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At B ^ Bt Eq. (|g) holds. Therefore, for c-axis uncor- 
related pancakes at T <S^ 70, aq{B,T) increases quasi- 
linearly with B with the slope 1/Bj when B <^ B^, and 
with a smaller slope ~ I/Sq, for B^ <^ B <^ Ba- At 
higher temperatures, T > 70, aq{B,T) decreases with 
field for B < Bt, before crossing over to the slow lin- 
ear growth with the slope ^ l/-Boj see Fig. 3. Here 
the conductivity has been computed with a temperature- 
independent Aq. In fact Aq is reduced with increasing T, 
and the plateau values increase. Note that linear interpo- 
lation of a{B, T) from high fields does not yield crq(0, T). 

To conclude, we find that: a) the field dependence 
of the quasiparticle interlayer conductivity is sensitive 
to the structure of vortex state; b) for the c-axis un- 
correlated vortices the field dependence of the interlayer 
conductivity at low and high temperatures is quite dif- 
ferent at low fields but becomes similar in high fields 
B > max[BT,-B7]- 
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FIG. 1. The distribution functions P{px 
a;-coniponent of supermomentum. Here p^ is measured 
in units h/ao, where ao = (2<I>o/\/3-B)^^^ is the lat- 
tice constant of a triangular vortex lattice. P{Px) has 
been calculated using Langevin dynamics simulations of a 
2D hquid at different reduced temperatures t = T/2-kEo. 
The corresponding absolute temperatures obtained assuming 
A„i, = 2QQnm(l - T^ /T^)'^''^ are given in brackets. The 
distribution function of the Doppler shifts P(e) is related to 
P{px) by 'P(e) = {2ao/hvF)P{px — 2aoe/TiVF)- Inset: tem- 
perature dependence of the parameter C\ in Eq. (6). 
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FIG. 2. Magnetic field dependence of the normalized quasi- 
particle conductivity cr5(_B,0)/cr5(0, 0) at low temperatures 
r <S 7o for c-axis uncorrelated vortex state calculated using 
functions P{px) shown Fig. 1. Here we do not show the weak 
linear field dependence at B/Bj ^ 1, see Eq. (6). 
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FIG. 3. Field dependence of the quasiparticle conductivity 
at temperatures T > O.570, computed for the pancake liq- 
uid state with P{px) obtained by simulations (Fig. 1) and 
70 = 0.16£o ~ 30 K. 



